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Abstract 

We study the mass spectra, production and decay properties of the lightest supersymmetric 
CP-even and CP-odd Higgs bosons in models with spontaneously broken R-parity (SBRP). We 
compare the resulting mass spectra with expectations of the Minimal Supersymmetric Standard 
Model (MSSM), stressing that the model obeys the upper bound on the lightest CP-even Higgs 
boson mass. We discuss how the presence of the additional scalar singlet states affects the Higgs 
production cross sections, both for the Bjorken process and the “associated production”. The main 
phenomenological novelty with respect to the MSSM comes from the fact that the spontaneous 
breaking of lepton number leads to the existence of the majoron, denoted J, which opens new 
decay channels for supersymmetric Higgs bosons. We find that the invisible decays of CP-even 
Higgses can be dominant, while those of the CP-odd bosons may also be sizeable. 
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I. INTRODUCTION 


Unveiling the mechanism of symmetry breaking and mass generation constitutes one of 
the main goals in the agenda of upcoming accelerators, like CERN’s Large Hadron Collider 
(LHC) and the International Linear Collider (ILC). Precision electroweak data currently hint 
that the mechanism responsible for electroweak symmetry breaking jlj involves a weakly- 
coupled Higgs sector, as predicted by supersymmetry. Supersymmetry also stabilizes the 
Higgs boson mass against quadratic divergences, thus accounting for the hierarchy between 
the electroweak and the Planck scales in a technically natural way. A very exciting pos¬ 
sibility is that the experimentally observed USB neutrino masses and mixings [J have 
a supersymmetric origin Jf|. The key requirement for this to be possible is that R-parity, 
defined as R p = (— \y B + L + 2S (with S, B , L denoting spin, baryon and lepton numbers, 
respectively) be violated. The simplest way to do this is through a bilinear term in the 
superpotential. 

The resulting model is interesting in two ways. First it provides the most economical 
description of R-parity violation as a “perturbation” to the Minimal Supersymmetric Stan¬ 
dard Model: it may be taken as the reference R-parity violation model, which we may 
call RMSSM 0]. The model offers a minimal low-scale mechanism to generate neutrino 
masses, that successfully accounts for the observed pattern of neutrino masses and mixing 
B, @] 1 . In contrast to the seesaw mechanism, it makes well defined predictions that will 
be tested at upcoming colliders LHC/ILC, namely, the decay branching ratios of the light¬ 
est supersymmetric particle are related to the neutrino mixing angles measured in neutrino 
experiments Jl3 |. 

On the other hand the model also provides the simplest effective description of theories 
where the breaking of R-parity occurs spontaneously, like that of the electroweak gauge 
symmetry itself, due to the existence of non-zero singlet sneutrino vacuum expectation 
values (vevs) H, Q HQ. A general feature of models where neutrino masses arise 
from low-scale spontaneous violation of ungauged lepton number is that the lightest CP- 
even supersymmetric Higgs boson will have an important decay channel into the singlet 
Goldstone boson (called majoron) associated to lepton number violation Q|: 


h^JJ. 


( 1 ) 


Thus the Higgs boson may decay mainly to an invisible mode characterized by missing 
energy, instead of the Standard Model channels. This general possibility can also be realized 
in spontaneously broken R-parity supersymmety Q]. 

We have recently reanalyzed this sug gest ion in view of the data on neutrino oscillations 
that indicate non-zero neutrino masses (l2|. We found that this proposal remains valid, 


1 For papers dealing with trilinear terms see, for example 0H 


2 







despite the smallness of neutrino masses required to fit current neutrino oscillation data [,5]. 
In Ref. 3 we have shown explicitly that the invisible decays of the lightest CP-even Higgs 
boson can be dominant, unsuppressed by the small neutrino masses, for the same parameter 
values for which Higgs production in e + e~ annihilation is comparable in cross section to 
that characterizing the standard case. A necessary ingredient in this case is the existence of 
an SU(2) <g) U(l) singlet superfield coupling to the electroweak doublet Higgses, which 
may provide a solution to the so-called //-problem. 

In this follow-up paper we extend the analysis and study in detail the possibility that 
the lightest CP-even Higgs boson is produced also in association with a CP-odd boson in 
electron-positron collisions. The first aspect to consider is the theoretically expected mass 
spectra of CP-even and CP-odd scalar bosons. An important feature of any supersymmetric 
model is the existence of an upper bound on the mass of the lightest CP-even scalar boson. 
We verify explicitly that this feature emerges in the present model. We also explain how the 
supersymmetric Higgs boson mass upper limit should be understood in terms of the SBRP 
model fields. 

Then we turn to the Higgs production cross sections. Although individual Higgs boson 
production cross sections, via the familiar Bjorken process or the associated mode, are 
potentially suppressed with respect to those of the MSSM, given enough center-of-mass 
energy and luminosity, all Higgses can be potentially explored due to unitarity. We carefully 
analyze the decay properties of the first and second lightest Higgs bosons. The case of the 
lightest CP-even Higgs boson was already considered in Ref. [12|. We revisit this case, and 
confirm that the invisible dacay mode for either the first or the second lightest CP-even 
Higgs boson can easily dominate, in contrast to that of the lightest CP-odd which may arise 
at subleading level up to 20 % level at most. 


II. THE MODEL 


For completeness we recall here the main ingredients of the model. In addition to the 
Minimal Supersymmetric Standard Model superfields it contains SU(2) ® U(l) singlet 


superfields (uf, S), <E>) carrying lepton number assigned as (—1,1,0). Wit 
most general superpotential terms conserving lepton number are given as 


l this choice the 



ti j $V?Sj + MffiSj + iM*8 2 + 


( 2 ) 


The first three terms together with the jl term define the R-parity conserving MSSM, the 
terms in the last row involve only couplings among the SU(2) ® U(l) singlet superfields 
(Pf, Si, <&). The remaining terms couple the singlets to the MSSM fields. We stress the 
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importance of the Dirac-Yukawa term which connects the right-handed neutrino superfields 
to the lepton doublet superfields, thus fixing lepton number. 

Like all other Yukawa couplings in general h u is an arbitrary non-symmetric complex 
matrix in generation space. However, for technical simplicity we will consider only the case 
with just one pair of lepton-number-carrying SU(2) ® U(l) singlet superfields, v c and S, 
in order to avoid inessential complication. This in turn implies, h u —> h and h l J —> h l u . 

The scalar potential along neutral directions is given by 


V total = \h<&S + hl,i>iH u + M r S\ 2 + \h 0 $H u + fiH u \ 2 + \h$v c + M R u c \ 2 

+ | - h 0 $H d - fiH d + + | - h 0 H u H d + hv c S - b 2 + + ^<f> 


( 3 ) 


212 


\K" c Hu\ 2 + - A ho h 0 $H u H d + A h X"iH u v c - BfiH u H d 


i =1 


1 „ - ' -2 I 3 A x 1 3 


-C s S 2 $ + B Mr Mrv c S + + -Ax + h.c. 
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+ > m a \z a \ 


+ g (g 2 + g' 2 ) ( \H U \ 2 - \H d \ 2 - 


i= 1 


where z a denotes any neutral scalar held in the theory. For simplicity we assume CP 
conservation in the scalar sector, taking all couplings real. 

Electroweak symmetry breaking is driven by the isodoublet vevs (H u ) = and (H d ) = 
with the combination v 2 = v 2 + v‘ d + Yli v u hxed by the W mass, while the ratio of 
isodoublet vevs yields tan/5 = Here, Y| are the vevs of the left-scalar neutrinos. They 
vanish in the limit where h l u —> 0. In this limit R-parity is restored and neutrinos become 
massless, as in the MSSM, and, apart from <f>, the extra singlets become phenomenologically 
irrelevant, one reaches the NMSSM limit 211. 

The spontaneous breaking of R-parity is driven by nonzero vevs for the right-scalar neu¬ 
trinos. The scale characterizing R-parity breaking is set by the isosinglet vevs (V c ) = ^ 
and (S) = ^|. Finally, ($) = gives a contribution to the /i-term. 

With the above choices and definitions we can obtain the neutral scalar boson mass 
matrices as described in Ref. [l4]. This results in 8 x 8 mass matrices for the real and 
imaginary parts of the neutral scalars. Their complete definition can be found in flo| . 
The spontaneous breaking of SU(2) 0 U(l) and lepton number leads to two Goldstone 
bosons, namely G° , the one “eaten” by the Z°, as well as J, the majoron. In the basis 
P'° = (H® 1 , H® 1 , v 11 , v 21 , z> 3/ , S 1 , v cI ) these helds are given as, 


G° — (N 0 v d ,—N 0 v u , N 0 vli, N 0 vl 2 , N 0 vlz,0,0,0) 

J = AT 4 (-iViu d , NxVu, N 2 v l i, N 2 v L2 , N 2 v l 3 , 0, N 3 v s , -N 3 v R ), (4) 
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where the normalization constants A) are given as 

No = 


V v d + v l + v'h + v 2 L2 + 


Ni — v L1 + v L2 + vl 3 
N 2 = V 2 d + vl 


N% — Ni + N 2 
N a = 


( 5 ) 


y/N*N 2 + iVfAh + N 2 (v 2 R + V 2 S ) 

and can easily be checked to be orthogonal, i. e. they satisfy G° ■ J = 0. 

The neutrino masses and mixings arising from this model [3] have been shown to repro¬ 
duce the current data on neutrino oscillations that indicate non-zero neutrino masses Js]. 
Since neutrino masses are so much smaller than all other fermion mass terms in the model, 
once can find the effective neutrino mass matrix in a seesaw-type approximation. After 
some algebraic manipulation, the effective neutrino mass matrix can be cast into a very 
simple form 

( 6 ) 


')*?' — uAjAj + b(etAj + e^-Aj) + ce 




where one can define the effective bilinear R-parity violating parameters e* and A* as 

Vr 


and 


Here the parameter /i is 


— h l 

1 U V2 

K = 6iV d + flV Li 


h — A + ' l 0 , 


( 7 ) 


( 8 ) 


( 9 ) 


while the coefficients appearing in Eq. © are given in Ref. jlflj | . Eq. © resembles closely the 
structure found for the explicit bilinear model at the one-loop level. However, the coefficients 
are different, see 

Neutrino physics puts a number of constraints on the parameters A, and e t . For the 
current paper, however, exact details are unimportant, the most essential constraint for the 
following discussion is that hi -C 1 is required. (See later discussion of left-sneutrino mixing. 
In the limit h l v — 0 left-sneutrinos do not mix at all with Higgses and singlets). 

The requirement that Vl l v can be used to fold a simple approximation formula for 
the majoron, given by 


J ~ 


-Vdl’l V u V 2 l V L1 V L2 V L3 Vs Vr 

Vv 2 ’ Vv 2 ’ V ’ V ’ V ’ ’ V ’ V ] ’ 


( 10 ) 


where V 2 = v 2 s + Vr. Thus, the majoron is essentially made up of the v c and S helds. This 
will be important later, when we discuss the decays of the Higgs bosons. 
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III. HIGGS SPECTRUM 


Let us first briefly discuss the spectrum of the scalar and pseudo-scalar sectors in the 
model. For detailed definitions we refer the reader to Ref. |19|. Since these mass matri¬ 
ces are too complicated for analytic diagonalization, we will solve the exact eigensystems 
numerically. However, before doing that, we discuss certain limits, where some simplifying 
approximations are made. This allows us to gain some insight into the nature of the spectra. 

In the SBRP model there are 8 neutral CP-even states S®. In the neutral CP-odd sector 
there are six massive states P,° (i — 1,..., 6), in addition to the majoron J, with rrij = 0, 
and the Goldstone G°. We introduce the convention, to be discussed below: 

{S°) T = (Sfco, S H 0 , Sj, S J± , Sz, S 0i ) (11) 

(P°) T = (P A o,P J± ,P*,P Pi ,J,G°). 


Note, that the ordering of these states is not by increasing mass, as we have defined P® 
(i = 1 , ..., 6) as the massive states. 

First we note that all entries in the sub-matrices which mix the left-sneutrinos to the 
doublet Higgses and the singlet states are proportional to h l v . In the region of parameters 
where the model accounts for the observed neutrino masses we must have that e* = h l u Vn/y/2 
is necessarily a small number and therefore h l u ■ Cl. Thus, left sneutrinos mix very little 
with the other (pseudo-)scalars, unless entries in the sneutrino sector are, by chance, highly 
degenerate with the ones in the other sectors. The real (imaginary) parts of these nearly- 
sneutrino states are denoted by S (PpJ in the definition given above. Barring fine-tuned 
situations, we conclude that mixing between Higgses and left sneutrinos will, in general, be 
small. 

Consider now the pseudoscalar sector, 



'Mhh 

Kl 

Mfts 

M p2 = 

M p2 - T 

MP- 

MP 2 


HL 

LL 

LS 






where is a symmetric 2x2 matrix, AP£~ and Mgg are symmetric 3x3 matrices, 

while and Mf? s are 2x3 matrices and finally is (a non-symmetric) 3x3 matrix. 

In this notation L denotes the sneutrinos and S the singlet fields. 

Neglecting terms proportional to h l u , can be written as 


Mhh = 


Vd 

n ri— 

V u 


(13) 
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where 2 


0 


A 1 

Bfi - S 2 h 0 + -/io^l + -hh 0 v R v s + — A ho h 0 v^ + —h 0 M^ . 


(14) 


Note the presence of ho-dependent terms in Eq. m- If there were no mixing between the 
doublet and singlet Higgses, Eq. m would yield the eigenvalues, 



(15) 


with the massless state identified as the Goldstone boson, G°, and the other state as the 
pseudo-scalar Higgs A 0 of the MSSM, with 


m 2 A o 


20 

sin 2 (3 ' 


(16) 


The state most closely resembling the MSSM H°, i.e. the state remaining in the spectrum 
when singlets are decoupled is called P 40 in Eq. ED • The sub-matrix M P g, on the other 
hand, in the limit h l v = 0 , can be written as, 


where, 


M p2 = 
lvl ss 


Mss u Mf Sl2 M * SSl3 

MsL -r“ -r 
M sL -r -rf, 


M, 


p2 


SS 11 


M, 

M, 


= V (C, + M t ) Ak _ (3^ + M») b* 

Z C/<J> ht 




— 2 — h (A h + M$) — -|——/to (AU 0 + -^$) 


\/2 . 


2 

4 + 4 


+ 25 A + Ah 0 v u v d - A h v R v s —h M R 

2 


SS12 — M®) v R j 


SS1 3 — M®) Vs ■ 




Here + At>$ / v^2 and 


1 1 \/2 

r = B Mr M r - S 2 li + -hXvl - -hh 0 v u v d + —h (A h + M$) . 


(17) 


(18) 

(19) 

( 20 ) 


( 21 ) 


2 We correct a misprint in 
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Eq. (IT71) has one zero eigenvalue, approximately identified with the majoron, J, and two 
non-zero eigenvalues. If M RRl2 , M RR -C MsSn + E then the eigenvalues of Eq. (H71) are 
approximately given by 


m 


2 _ 

1 , 2,3 — 


n r( VR \ Vs ) 1 h 2 (A h - MzYv%v 2 s 
v s v R 2 Mf Sii v R v s + T(v 2 r + v 2 s ) 

m p 2 , 1 h?(A h - M^) 2 v r v s {v r + v 2 s ) \ 

SSl1 2 Mf Sll v R vs + r(4 + v 2 ) )’ 


( 22 ) 


where the dots stand for higher order terms. The eigenvalue proportional to T is mainly a 
combination of S 1 , D cI fields and we call it Pj ± in Eq. (1111) above, because in the limit where 
mp j, —> oo and v Ri —■► 0 this massive state is orthogonal to the majoron. As we will discuss 
below, it is this state which preferably decays invisibly. The third eigenvalue in Eq. E2D is 
an approximation to the state called above. Due to mixing between doublet and singlet 
states both Eq. m and Eq. E2D, are only very crude estimates. 

Consider the scalar sector of the model, 


M b = 


M HH Kl M hs 
mS h 7 Ks 


M 


S 2 T 


M i7 M ss 


rS 2 


HS 


(23) 


where the blocks have the same structure as before. M$h contains two eigenvalues which, 
in the limit of zero mixing, would be identified with the MSSM states h° and H°. 3 These 
states are the ones called S^o and S R o in Eq. m above. 

The sub-matrix contains, in general, three non-zero eigenvalues. One can find an 

approximate analytic expression for them in the limit that the state <S$ is much heavier than 
the remaining two eigenstates (called Sj and Sj ± ). Again in the limit of small mixing, the 
eigenvalues of the latter are approximately given by 


m 


1,2 


= 2 h 


V R V S 


J R 


+ «I) 


+ --,-T( 


VR 

Vs 


+ —) — 2h 2 

VR 


v 2 R v S 


J R 


+ v 2 s) 


+ ••• 


(24) 


The first (second) of the eigenvalues in Eq. EH) is approximately the state Sj ( Sj ± ). 

Fig- CJ to the left, shows an example of the four lowest lying eigenvalues in the CP-even 
sector, as a function of T for a random but fixed choice of the remaining parameters. One 
of the states, Sj ± , which is mainly singlet, is proportional to T, as indicated by Eq. EH) - 
There is another singlet state, corresponding to Sj of Eq. EHb an d two mainly doublet 
states, identified with S h o and S H o. We note in passing that rns H0 i s proportional to O, as 
in the MSSM. Mixing between singlet and doublet states will be important always if the 


3 As in the MSSM, there is an upper limit for the mass of the Sw , see the discussion below. 









eigenvalues are comparable, as for the example shown in the figure. Thus, all the discussion 
above should be taken as qualitative only. 

The right panel in Fig. Q] shows an example of the two lightest massive CP-odd eigenvalues 
as a function of T for a fixed but random set of other parameters. That one eigenvalue is 
proportional to T is obvious from Eq. (1221) . We note that and T are the main parameters 
which will determine associated production and influence the branching ratio into invisible 
states, as we will discuss in the following sections. The model clearly exhibits decoupling, 



Figure 1: Typical CP-even (left) and CP-odd (right) Higgs masses as function of the parameter 
r. In this example there are four light CP-even states and two light massive CP-odd states (plus 
two massless states, G° and J, not shown). Just as in the MSSM there is always one light doublet 
state, coinciding with h° in the limit of zero mixing. Other states can (but need not) be light, 
depending on the parameters H and T, see text. 

just as the MSSM. In the limit where O goes to infinity the masses of both states Pa° and 
Sh o go to infinity, just as what happens in the MSSM when m^o goes to infinity. The states 
Sj and Pj ± are decoupled in the limit as V goes to infinity. If, in addition, we require 
h <C 1 also Sj decouples and the SM Higgs phenomenology is recovered, as in the MSSM. 


IV. HIGGS BOSON PRODUCTION 

Supersymmetric Higgs bosons can be produced at an e + e _ collider through their couplings 
to Z °, via the so-called Bjorken process (e + e~ —> Z°Sf), or via the associated production 
mechanism (e + e _ —> SfPj). In our SBRP model there are 8 neutral CP-even states Sf and 
6 massive neutral CP-odd Higgs bosons Pf, in addition to the majoron J and the Goldstone 
G°, see Eq. (fill) . 

One must diagonalize the (pseudo-)scalar boson mass matrices in order to find the cou- 
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plings of the scalars to the Z°. After doing that we obtain the Lagrangian terms 


C D £(V 2G f ) 1/2 MlZjZ°^ Bi S? + ^(V2 G F ) 1 / 2 M zVAll (z a >‘S'?%P?) (25) 

2—1 Z,J = 1 

with each t/b, given as a weighted combination of the five SU(2) ® U(l) doublet scalars, 


and the r)\ i] given by 


„ _ pS° , V u r,S° 

lBi ~ u 11 + ^ i2 


v Ri i+ 2 

3 = 1 


(26) 


r, A „ = + f? R? m R? m (27) 

fc=l 

where the subscripts B and A refer, respectively, to the Bjorken process or associated pro¬ 
duction mechanisms. From these Lagrangian terms we can easily derive the production cross 
sections. These are simple generalizations of the MSSM results 122i. |23|| and for completeness 
we give them in Appendix 0 

In the MSSM, there are two sum rule rules, one concerning only the CP even sector 



+ vl 


H 0 


= 1 , 


and another relating the Bjorken and the associated production mechanisms, 


(28) 



+ ’)a 


hOAO 


1 , 


(29) 


with T]B h o = sin(a — (3 ) and Va . 0a0 = Vb h0 — cos(cx — / 3 ), in an obvious notation. 

It is very easy, and instructive, to use our expressions for rj a and r )b to recover the MSSM 
result in the limit that we have only the H^ and H u doublets. In fact, in this case 


Vd 

v 


R 


pO 
22 ) 


Vu 

V 


R 


P° 

21 


so we have 


(30) 


„ _ nP° tdS° , pP° r S° _ pP° p5° I pP° p5° 

VB h 0 — Vl‘2'2 -till I 0-21 rt 12 , HB h o — ■'*'22 'Ml +'Ml -''22 ) 


= <R S u 


R%Rf 2 ( 31 ) 


and we get for the sum rule of Eq. (ESI) 


vl 


hOA® 


+ vl 


hO 


l R p0 R-S 10 I pP° R s°\ 2 I ( tdP° jdS° _ r>P° p S°\ 2 

I Jl 2 2 -tin \ ±i21 ^12 ) ‘ l "21 -^11 -^22 "12 J 

(a?; a?; + R?? a??) (it?? it?? + rSrS) 

i, 


(32) 

(33) 

(34) 
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where we have used the orthogonality of the rotation matrices 


= '£, R * R % = 6 <i’ (>,i = 1,2). 


fc=l 


k=1 


For the sum rule of the CP-even sector, Eq. ( 121 , we get 


’Jb.o+’Jb.o = p [ Rf, Rf, + Fti; Ri;)+ sin 1 13 (Rf°Rf 2 +Rl°Fti° 


+2 sin/3 cos/J (flf°-Rf 2 " + R 2 ° R S 


= 1 


( 35 ) 


(36) 

(37) 

(38) 


using the result that in an orthogonal matrix also the vectors corresponding to the columns 
are orthonormal, that is 

2 

J2 R l‘ R l‘ = S ‘i’ (Uj = h 2). (39) 

fc=l 

How this differs in our case? The difference is that, in general, 

2 2 

E T}S°,P° tjS°,P° / r J X"'' r,5°,P° T}S °, p0 _L c liol 

/c=l /c=l 

due to the fact that we now have more than two (pseudo-)scalars. As it was stated in the 
last section and will be discussed in more detail when we consider the decays, to have a 
sizeable invisible branching ratio we need the doublets to be close in mass to the singlet 
states related to the majoron and orthogonal combinations. This means that, in the CP- 
even sector, the first four states are (5\o, Sh°, Sj ± , Sj ), while in the CP-odd sector we should 
have (-Pa 0 , A G°). If this situation happens then we can very easily find a generalization 
of the sum rule of the CP-even sector, as 


2 , 2 , 2 , 2 
+ 7 1b S h0 + Vb S j± + Vb s 


= 1 


(41) 


to a good approximation. This is displayed in Fig.Owhere we plot the sum rfy +?7jL + 

S H 0 S J± 

against . The significance of this sum rule should be clear: if the lightest Higgs 
boson has a very small coupling to the Z° and hence a small production cross section, there 
should be another state nearby that has a large production cross section. 

The other sum rule, relating the CP-even and CP-odd sectors, Eq. (EHD, is more difficult 
to generalize. In fact the Pa 0 state will now mix with the Pj ± and the identification of 
Eq. (T¥)l) will be no longer true. However qualitatively the sum rule still holds in the sense 
that if the parameters are such that the production of the CP-odd states is reduced one 
always gets a CP-even state produced. 
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Vbs 


h° 


Figure 2: Sum rule in the CP-even sector, for the case explained in the text. The four states, 
(5/jO, S B o, Sj ± , Sj). For this example all scalar masses are taken below 200 GeV. 


The above discussion has concentrated on Higgs boson production at an e + e collider. 


We now briefly comment on the differences with regards to Higgs production at the LHC 24 ]. 


It has 
fusion 


aeen suggested to search for an invisibly decaying Higgs at the LHC in WW boson 


rr 


in asociated production with a Z° boson j26(], or in the ti channel |27[. For the 
production in WW fusion or in asociated production with a Z° boson the above discussion 
applies straightforwardly, since the relevant coupling in both cases is r] B . (i.e. sin{j3 — a) in 
the MSSM limit). For the ti channel in the MSSM production cross section the factor cos a 
has to be replaced by for the SBRP model. 


V. HIGGS BOSON DECAYS 


In the following we will discuss the decays of light CP-even and CP-odd supersymmetric 
Higgs bosons. Since the phenomenology of Higgs bosons within the MSSM is well-known 0 
1291 ]. we will concentrate on non-standard final states. Of these, the most important are the 
majoron Higgs boson decay modes, which are characteristic of the SBRP model, without 
an MSSM counterpart. We will limit ourselves to the discussion of light states, i.e. Higgs 
bosons with masses below the 2W threshold. As discussed below, the decays of heavier 
CP-odd states will be similar to the situation encountered in the (N)MSSM. 


A. CP-even Higgs Boson Decays 

In the MSSM light CP-even Higgs bosons decay dominantly to bb final states. In our 
calculation we take into account all fermion final states, including the leading QCD radiative 
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corrections from 


if kinematically allowed S? 


30]. In the SBRP model new decay modes appear, such as Sf —> JJ and, 

P?Pu- From the latter usually only S 0 —> JJ 


poj and Sf 

J L 


has a large branching ratio (see appendix El) • 

In Ref. fig] we have discussed the invisible decays of the lightest CP-even Higgs boson. 
We now extend that discussion so as to include also the next-to-lightest CP-even state which 
plays an important role, if the lightest CP-even state is mainly singlet. 

It is well known that, in contrast to the Standard Model, in the MSSM (and in the 
NMSSM) the mass of the lightest CP-even supersymmetric Higgs boson obeys an upper 
bound that follows from the D-term origin of the quartic terms in the scalar potential, 
contained in Eq. (l3l) . This mass acquires a contribution from the top-stop quark ex¬ 


change 


31 


32, 


, a fact that modifies the numerical value of this upper bound 


31 


0 


0 , 0 . 

This 


Many other loops contribute, for a recent two-loop level calculation see, e.g. Ref. 
limit is slightly relaxed in the NMSSM as opposed to the MSSM |0 . 

How does this bound emerge in the SBRP model? Since the CP-even sector contains eight 
scalars, we cannot diagonalize the corresponding mass matrices analytically. Therefore we 
calculate the upper bound on the Higgs mass numerically, and including the most important 
radiative corrections, using formulas from 0 . In the SBRP model it is possible that the 
lightest CP-even Higgs is mainly a singlet. However, if this happens, there must exist a 
light, mainly doublet Higgs, to which the NMSSM bounds apply. This is shown in Fig. 01 
where we plot (to the left) r /^ 2 as function of the 77 ^ and (to the right) the upper limit on 
the mass of the second lightest Higgs as function of ?/| 0 . As is seen, if the lightest state 
is mainly singlet, 77| ~ 0 , therefore r /^ 2 ~ 1, then there is an upper bound on the second 
lightest state mass. Vice versa the upper bound applies to the lightest state if it is mainly 
doublet. 

As shown previously [0, one can have large direct production cross section for the 
lightest neutral CP-even Higgs boson as well as a large branching ratio to the invisible final 
majoron states. This is demonstrated in the left panel of Fig. E for a random but fixed 
choice of undisplayed parameters. We note that a very similar behaviour is also found for 
the second lightest state, as seen from the right panel of Fig. [U Thus if the lightest state is 
mainly singlet there must be a state nearby which is mainly doublet and decays invisibly. 

In summary, we have seen that in the SBRP model there is always at least one light state, 
which is mainly doublet, and therefore can be produced at future colliders. Irrespectively of 
whether this state is the lightest or second-lightest Higgs state, it can decay with very large 
branching ratio to an invisible final state. 
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Figure 3: In the left panel we show the parameter characterizing direct production of the second 
lightest neutral CP-even Higgs boson, , as function of the corresponding one for the first lightest 
neutral CP-even Higgs boson, . To the right: Upper limit on the mass of the second lightest 
CP-even Higgs as a function of 7?| o . 
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Figure 4: To the left (right): Ratio R± (R 2 ) as a function of the direct production parameter, r/jy 
(? 7 g 2 ), for the first (second) lightest neutral CP-even Higgs boson. 


B. CP-odd Higgs Boson Decays 

Light CP-odd Higgs bosons in the MSSM decay according to Pf —> //. The WW 
channel becomes dominant as soon as kinematically allowed CH Eoj]. however we will not 
include it as we are mainly interested in the possibility of invisible decays of the lowest- 
lying pseudoscalar. The formulas for the CP-even and CP-odd Higgs boson MSSM decay 
branching ratios, apart from the larger number of Higgs bosons, are totally analogous to 
those of the MSSM j^3|, except for the prefactors which are determined by the diagonalizing 
matrices of our model. The corresponding matrix elements replace the familiar sin (/3 — a) 
and cos (/3 — a) factors. 
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In the SBRP we must take into account in addition the decays P} } —> JJJ and, if 
kinematically allowed, also P f —■> SjJ, P t ° SjP%, P ? ° —■> PfJJ , Pf — > PjP%J and 

Pf —> PjP®P® 2 . For the lightest Higgs boson we are interested only in P° —> JJJ and 
P° —> SjJ. The formulas for the CP-even and CP-odd Higgs boson non-MSSM decay 
widths are collected in appendix |Bj 


C. P° 


S«J 


The decay width of the CP-odd Higgs boson to a CP-even Higgs boson and a majoron 
is given in Eq. (Hu. Using the approximation Eq. (HDD we can find the coupling g' %rj for the 
vertex S'^P'^J of the majoron with the unrotated neutral scalar S® and pseudoscalar Pj° to 
leading order in the small parameter — as 
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(42) 


Note that the first four of the above are suppressed by the smallness of sneutrino vevs, needed 
to reproduce the observed neutrino oscillation data. The coupling gs°p 9 j then appears 

1 3 
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through mixing, and is given as 


asfpfj = at« + + a' m R?sRn + 9^« 2 ' 


S° D P° 


?S° D P° 


?S° dP° 


(43) 


D. P° JJJ 


The decay width of the CP-odd Higgs boson to three majorons is given in Eq. (ESI). Using 
again the approximate equation giving the profile of the majoron, Eq. (DUD, the coupling g[ 
for the vertex P-°JJJ of the majorons with the unrotated neutral pseudoscalar P/°, is given 


as 
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(44) 


Again, the first six of the above vanish in the limit v R —> 0. Therefore the coupling gpajjj 
for the vertex of the majorons with the neutral pseudoscalar P )° mass eigenstate is 


_ / D p° . / D p° 

9 p?jjj — 97^17 + 9 %Pis ■ 


(45) 


E. Numerical results 

We can see from Eq. (H21l that if the CP-odd mass eigenstate is mainly a Higgs doublet 
(i.e., its main components are P[° = H® 1 , P = H°J so that its production is not reduced) 
then its decays to SjJ and JJJ are suppressed as the corresponding couplings are very 
small, suppressed by two powers in —. To find sizeable branching ratios for the decays of 
the lightest massive pseudoscalar P,°, mixing between doublet and singlet states is therefore 
required. 

As discussed in section nn in order to have sizeable mixing between doublet and singlet 
CP-odd Higgs bosons, we must require that at least one of the singlet states is light, i.e. 
the parameter V should be very roughly of order V rv./ fl Fig. El shows an example. Here, 
we plot r]\ 21 and BR(P 1 ° —> inv) as function of for one fixed, but arbitrary set of other 
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Figure 5: Production cross section (red/solid curve) and invisible final states decay branching ratio 
(green/dashed curve) for the lightest CP-odd Higgs boson. 

model parameters. For small values of y/T the lightest massive CP-odd state is mainly 
singlet, therefore BR(P 1 ° —> inv) is close to 1. However, the production parameter r}\ is 
small. Increasing \/r increases the mass of the lightest CP-odd state. From a certain point 
onwards, it is the doublet state which is lightest, compare to Fig. Q] This state can have a 
sizeable production, but the branching ratio to invisible final states typically is small. Only 
in the intermediate region of sizeable mixing between doublet and singlet states, i.e. in the 
region of \/r rs j 100 — 115 GeV of Fig. 0 one can have both, sizeable production and sizeable 
invisible decay. 

In summary, the CP-odd Higgs bosons in the SBRP model usually behave very similar 
to the situation discussed in the (N)MSSM. However, sizeable branching ratios to invisible 
final states are possible when there are light CP-odd Higgs bosons from both, the doublet 
and the singlet sectors. 


VI. DISCUSSION 

We have carefully analyzed the mass spectra, production and decay properties of the 
lightest supersymmetric CP-even and CP-odd Higgs bosons in models with spontaneously 
broken R-parity. We have compared the resulting mass spectra with what is predicted in 
the Minimal Supersymmetric Standard Model, stressing the validity of the upper bound on 
the lightest CP-even Higgs boson mass. We have seen how the presence of the additional 
scalar singlet states affects the Higgs production cross sections, both in the Bjorken and 
associated modes. 

The main difference with respect to the MSSM case comes from the fact that the sponta¬ 
neous breaking of lepton number necessarily implies the existence of the majoron, and this 
opens new decay channels for supersymmetric Higgs bosons into “invisible” final states. We 
have found that the invisible decays of CP-even Higgses can be dominant, despite the small 
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values of the neutrino masses indicated by neutrino oscillation data. In contrast, although 
the decays of the CP-odd bosons into invisible final states can also be sizeable, this situation 
is not generic. 

Therefore the existence of invisibly decaying Higgs bosons should be taken seriosly in 
the planning of future accelerators, like the LHC and the ILC. These decays may signal the 
weak-scale violation of lepton number in a wide class of theories. Within the supersymmetric 
context they are a characteristic feature of the SBRP models. These can account for the 
observed pattern of neutrino masses and mixings in a way which allows the neutrino mixing 
angles to be cross checked at high energy accelerators like LHC/ILC. For example, in our 
model there is a bb plus missing momentum signal associated to the invisible decay of the 
lightest CP-even Higgs boson produced in association with a pseudoscalar. Although this 
is a standard topology, also present in the Standard Model and the MSSM, its kinematical 
properties in our model differ, as the JJ add up to the CP-even Higgs boson mass and bb to 
the CP-odd Higgs boson mass. Further studies to elucidate the impact of these decay modes 
for future colliders, should be conducted. While for the LHC we may encounter difficulties 
associated to missing energy measurements and/or b-tagging, these potential limitations do 
not affect in the same way the ILC. 

Last, but not least, as already explained, we have restricted our analysis to Higgs bosons 
below the WW threshold. Extension to relax this restriction is totally straightforward, 
though somewhat less interesting. Due to the validity of the supersymmetric Higgs boson 
mass upper limit we must have one light CP-even Higgs boson which, as we have shown, is 
likely to have an important decay into invisible final states. 
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Appendix A: PRODUCTION CROSS SECTIONS 

In this section we give the formulas for the production cross section of both channels at 
an e + e“ machine. 
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1. Bjorken process 


The cross section for the Bjorken process is 221 

G 2 M% 


a(e + e -> Z°S °) = r/^ {v 2 e + a 2 e ) (3 


/3 2 + 12M z /s 


where 


v e = -1 + 4 sin a e = -1, (3 = 


(1 - M 2 /s) 2 + {YzMz/sf ’ 
A(s,Mf,Mj 0 ) 


A is the 2-body phase space function, 


A (a,b,c) = yj(a + b — c) 2 — 4 ab 
and the ?/b, are given in Eq. (126(1 . 


(Al) 

(A2) 


(A3) 


2. Associated production 


The cross section for the associated production is 
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and the rj are given in Eq. (EZD- 
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Appendix B: NON-MSSM DECAYS 

The most characteristic decays of this model which do not exist in the (N)MSSM are 
those involving a majoron. In the following we collect the formulas for these decays. The 
most important ones are: 
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For completeness we consider also: 
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The decays Pj° —> P°P°P° are possible, but closed kinematically for the light states of 
interest, therefore we do not give here the explicit formulas for the widths. 
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